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Abstract 

We have studied the evolution of the massless scalar field propagating in time-dependent charged Vaidya 
black hole background. A generalized tortoise coordinate transformation were used to study the evolution 
of the massless scalar field. It is shown that, for the slowest damped quasinormal modes, the approximate 
formulae in stationary Reissner-Nordstrom black hole turn out to be a reasonable prescription, showing that 
results from quasinormal mode analysis are rather robust. 
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I. INTRODUCTION 



It is well known that the surrounding geometry of a black hole will experience damped oscil- 
lations under perturbations. These oscillations are called "quasinormal modes" (QNM), which is 
believed to be a characteristic "sound" of black holes, and would lead to the direct identification of 
the black hole existence through gravitational wave observation, to be realized in the near future. 
In the past few decades, a great deal of effort has been devoted to the study of QNMs of black 
holes in asymptotically flat spacetimes (for comprehensive reviews see [1,2] and references therein). 
Considering the case when the black hole is immersed in an expanding universe, the QNMs of 
black holes in de Sitter space have also been investigated [3-5]. Motivated by the recent discovery 
of the Anti-de Sitter/ Conformal Field Theory (AdS/CFT) correspondence, many authors have 
performed the study of QNMs in AdS spaces [6-18]. The study of the QNMs plays an important 
role in astrophysics, black hole physics and also string theory. 

All these previous works on QNMs have so far been restricted to time-independent black hole 
backgrounds. It should be realized that, for a realistic model, the black hole parameters change 
with time. A black hole gaining or losing mass via absorption (merging) or evaporation is a good 
example. The more intriguing investigation of the black hole QNM calls for a systematic analysis 
of time-dependent spacetimes. Recently the late time tails under the influence of a time-dependent 
scattering potential has been explored in [19], where the tail structure was found to be modified due 
to the temporal dependence of the potential. The motivation of our work is to explore the modifi- 
cation to the QNM in time-dependent spacetimes. Instead of plotting an effective time-dependent 
scattering potential by hand as done in [19], we will introduce the time-dependent potential in a 
natural way by considering dynamical black holes, with black hole parameters changing with time 
due to absorption and evaporation processes. We will study the temporal evolution of massless 
scalar field perturbation. Instead of employing Kruskal-like coordinate in our first attempt before 
[20], here we start our discussion directly from the Vaidya metric. 

The outline of this paper is as follows. In Sec. II we first go over the conventional treatment for 
the study of the wave propagation in stationary Reissner-Nordstrom (RN) black hole. In Sec. Ill, 
we will present the numerical way for the study of the QNMs of charged Vaidya black hole. We 
will present our numerical results in Sec. IV. Sec. V contains a brief summary and a discussion. 
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II. QUASINORNAL MODES OF A STATIONARY REISSNER-NORDSTROM BLACK 
HOLE 

Perturbation outside the stationary RN black hole has been discussed by many authors [21]. 
Here we first present a brief review on the conventional treatment of the evolution of massless 
scalar field propagating in RN black hole background. 

The metric of a RN black hole in the ingoing Eddington coordinates is given by 



where v is a coordinate usually called the 'advanced time', which parameterizes a lapse of time. 
The propagation of scalar waves in curved spacetimes is governed by the Klein-Gordon equation 



Since the background is spherically symmetric, each multipole of the perturbing field evolves sep- 
arately. Hence one can define 



ds 2 = -(1 - 2M/r + Q 2 /r 2 )dv 2 + 2dvdr + r 2 (d6 2 + sin 2 6d(p 2 ) = g^dx^dx' 



(1) 




(2) 




(3) 



Using the tortoise coordinate r* defined by 




(4) 



r + — r. 



r + — r 



where, r + = M + \J M 2 - Q 2 and r_ = M - \J M 2 - Q 2 are the radius of outer and inner horizons 
of the black hole, the wave equation for each multipole moment becomes 



# r „ r , + 2* )r ,„ - V* = 0, 



(5) 



where 





Using the null coordinate 



u = v — 2r* 



(7) 



equation (JHJ can be written as 



(8) 
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The two-dimensional wave equation (JSJ can be integrated numerically, using for example the finite 
difference method suggested in [22]. Using Taylor's theorem, it is discretized as 



N 



+ - *s - 6u5vV( 



VN + V W ~ U N - UE, + 



0{e% 



0) 



where the points N, S, E and W from a null rectangle with relative positions as: N : (u + 5u, v + 
5v),W : (u + 5u,v),E : (u,v + 5v) and S 1 : (u,v). The parameter e is an overall grid scalar factor, 
so that 5u ~ Sv ~ e. For an RN black hole in de Sitter space this has been performed in [4]. 

In addition to solving @ directly by numerical method, other ways have been suggested in 
studying the QNMs of RN black hole. For the slowest damped QNMs, the WKB formulas suggest 
the approximate behavior in RN black hole [23] 
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in the limit I » 1. Here ro = §M + \{9M 2 —&Q 2 ) 1 / 2 is the position where the potential V attains 
its maximum value. Defining q = Q/M, the above equation becomes 
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The WKB computation has also been performed in the RN de Sitter Spacetime [4]. 

These conventional treatments are powerful in studying the wave propagation in stationary 
black hole background, however it is difficult to extend it to the time-dependent case. For the 
charged Vaidya black hole, the mass and the charge are functions of time, M = M(v),Q = Q(v), 
and the standard tortoise coordinate can not be used to simplify the wave equation. 



III. FIELD EVOLUTION IN CHARGED VAIDYA BACKGROUND 

The metric of the charged Vaidya spacetime reads 

dg2 = _ 2M(u) + Q^l )dv 2 + 2cdydr + r 2 {dQ 2 + gin 2 Qd 2)) (12) 

where M = M(v) and Q = Q(v) are arbitrary functions of the time. For c = 1 the field is ingoing 
and M is monotone increasing in v (advanced time), and for c = — 1, the field is outgoing and M 
is monotone decreasing in v (retarded time) [24]. The QNMs of a black hole correspond to the 
solutions of the Klein-Gordon equation that satisfy the causal condition that no information could 
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leak out through the event horizon of the black hole and at the same time correspond to purely 
outgoing waves at spatial infinity. For the charged Vaidya black hole, horizons r± can be inferred 
from the null hypersurface condition 'f tfl f )V = and f(r±,v) = 0. r±(v) satisfies the equation 

r\ - 2r±M + Q 2 - 2cr 2 ± r± = 0, (13) 

where r± = dr±/dv. 

Since the charged Vaidya metric is spherically symmetric, a generalized tortoise coordinate 
transformation can be introduced as 

■ ln(r — r + ) — ln(r — r_), = v — vq, (14) 



2k + x ~ r ' 2k_ 

which is similar to equation ©, and which we take as a simplifying Ansatz for our numerical 
equations upon an appropriated choice of the constant s k± , as well as the functions r±(v), which 



will be interpreted as the event and Cauchy horizons [28(, r+ = M + y M 2 — Q 2 + 2cr 2 ^_r + and 
r_ = M — J M 2 — Q 2 + 2cr 2 _f_, respectively. We take k+ and fc_ as adjustable parameters to be 
defined. The parameter vq is an arbitrary constant assumed to be zero now. From formulas ()14|) 
and (112)) . the Klein-Gordon equation for each multipole moment becomes 



(1 + e 2 )^,r»r* + 2c¥ r .„, + ei*, r , - VV = 0, (15) 

where 
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with 
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For the convenience of numerical calculation, the adjustable parameters £;+ and can be 
selected in such a way that the Klein-Gordon equation for each multipole moment becomes the 



standard wave equation near the horizons r + (vo) and r_(t>o)- We define k + and by the algebraic 
equations 



r--M-2cr_r_ 1 r+-r_ 1 _ -, , „ . 

^ fc_ C r+-r_ fc+ _ l-hZCr_ 
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(17) 



at v = vq. We impose for simplicity, that 



lim E2 = lim e% = lim V = 0. (18) 
r -> r±(«o) r^r ± (u ) r -> r±(v ) 



f — > v o v — > v — > 



For the static black hole with r + = r_ = (or M and Q are constants), we have the usual result 



r 2 

= r _ r _ , when (JUJ) boils down to © and (|T|)J) to ©. 



Analogous to the coordinate used in @, we make the variable transformation 

u = it(r*, u*), v = (19) 



where, the curve ^(r*,^*) =constant is determined by the equation 

dr* I+E2 



(20) 



dv* 2c 

For £2^0 and c = 1, we have u — > i>* — 27"*, which is similar to l(7j). Using the transformation 
(|T^|). equation (|15|) can be simplified to 

5> - (2 C ^)- 1 F^ = 0. (21) 

The function u(r*,v*) can be integrated numerically according to equation (1201) . The integration 
of Vl/ can proceed similarly to equation l|9jl. 

However, the numerical calculation discussed above has some difficulty for small Q, especially for 
the Vaidya metric with Q — > 0. In this case the inner horizon r_ disappears and the transformation 
Q14|) is no longer valid. In this case (|14|> can be replaced by [25] 

1 

r * = r + TT7 m ( r ~ r +)> v * = v ~ v (22) 
2k 

where k is an adjustable parameter. Using (|22|). the Klein-Gordon equation for each multipole 
moment has the same form as equation (jl5j) . only replacing the expressions of £2, £i and V by 

£ 2 = [2k(Q 2 - 2Mr) + r + r + - 2M - 2c(r + r + )r + ]/{2kr 2 ), 
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The parameter A; can be selected to simplify the Klein-Gordon equation. We choose it such that 
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In the numerical calculation for the field evolution ^f, the coordinate transformation (122fi can be 
appropriated for Q — > 0. 

In the following numerical calculation for the field iff, we always use a 5 pulse as an initial 
perturbation, located at vo = 0, r = 3. We first use our numerical method to the case of stationary 
RN black hole. The numerical result is in agreement with equation (|1U|) or (jllj) . even for 1=2. 
Table 1 shows the slowest damped quasinormal frequencies for the scalar perturbation in the RN 
background with the multipole index 1 = 2. The numerical result is also in agreement with Ref. 
[4] for q=0. 

Table I: The numerical result for Mujr (or Mu>i) and the lowest WKB approximate formulas for cr(2, q) (or 
ci(2,q)) in the slowest damped quasinormal frequencies in the RN background. The multipole index / = 2. 
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c R (2,q) 
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Mujr [4] 


Mui [4] 





0.483 


0.481 


-0.0965 


-0.0962 


0.484 


-0.0965 


0.7 


0.532 


0.530 


-0.0985 


-0.0981 






0.999 


0.626 


0.624 


-0.0889 


-0.0886 







IV. NUMERICAL RESULT OF QNMS IN THE CHARGED VAIDYA BLACK HOLE 

It is widely believed today that asymptotically flat black holes eventually evaporate due to 
the emission of Hawking radiation. For a RN black hole, we assume that at a certain retarded 
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moment v = vq a charged nullfluid of negative-energy-density starts falling into the black hole 
(corresponding to the beginning of the evaporation process). We first adopt the 'linear model' [26], 
where M and Q depend linearly on v, and we also assume that the charge to mass ratio remains 
fixed. 



M(v) 



niQ v < vq 

mn(l — Xv) vn < v < V\ 

(26) 

mi v > v\ 



Q(v) = qM(v), 

where mo, mi, A and q are constant parameters. The above model is initially a RN black hole for 
v < vq. The evaporation process is displayed in the region vq < v < v\. At v = vx, the evaporation 
process stops, and the final geometry v > V\ is again the RN black hole. Changing the sign before 
A in (|26|) . this model can mimic the black hole absorption process. In the numerical study, we only 
consider how does the absorption or evaporation process affect the QNMs of the black hole in the 
region vq < v < V\. 

For the linear model (|2f?)) , we can get the solution r± of equation (fT3|) as 



r + = M + My/l - q 2 v>vi, 

M+MJl-q 2 (l-2r+) . 

T + = l-2r+ V ^ U l' 

M-MJl-q 2 (l-2r + ) 

r - = 1^2U « V ^ U 0, 

r_ = M — My/1 - q 2 v < v , 
for absorption process with c = 1, and 



(27) 



M+MJl-q 2 (l+2r + ) 
r + = , " > VQ, 



r+ = M + Ma/1 - q 2 v < v , 
r- = M - M v 7 ! - q 2 v > v 1} 



(28) 



M-MJl~q 2 (l+2r+) . 

r - = T+27T ' V ~ Ul ' 

for evaporation process with c = — 1. In general we suppose that 
to r± w M ± Ma/1 - g 2 . Thus the r + is very close to the outer apparent horizon M + Ma/1 — q 2 



M 



Amo << 1, which leads 



[26], while the location of the r_ is very close to the inner apparent horizon M — M\/l — q 2 . If 
q — > 1, then r + ~ M, r + ~ M. Meanwhile equations Q27|) and (|28|1 require l — g 2 (l — 2cr+) > 0, 
which is obviously satisfied for both the absorption and evaporation cases. As an example, we take 
mo = 0.5, A = 0.002, vo = 0, v\ = 150 (correspondingly, mi = 0.35 and M = —0.001). In Fig.l we 
show behaviors of r + and r_ for three different q in the evaporation process. 
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Figure 1: The behavior of r + and r_ for the black hole with q=0, q=0.7 and g=0.999, respectively. The linear 
model is used with mo = 0.5, A = 0.002, w = 0, uj. = 150 (correspondingly, mi = 0.35 and M = —0.001). 

The model ()26(l can represent the absorption process of the black hole when A < (or mo < mi). 
As an example, the temporal evolution of the field in the Vaidya metric (q=0) at r = 5 is displayed 
in Fig. 2. The relevant parameters were chosen as mo = 0.5, A = ±0.002, vo = 0,v\ = 150, I = 2. 
The initial 5 pulse perturbation is located at vq = 0, r = 3. For comparison, we also exhibit the 
curve obtained for Schwarzschild black hole (replacing the time-dependent mass by a constant mass 
M = 0.5). The modification to the QNMs due to the time-dependent background is clear. When 
M increases linearly with v, the decay becomes slower compared to the stationary case, which 
corresponds to saying that \u>i\ decreases with respect to v. The oscillation period is no longer 
a constant as in the stationary Schwarzschild black hole. It becomes longer with the increase of 
time. In other words, the real part of the quasinormal frequency lor decreases with the increase 
of time. When M decreases linearly with respect to v, compared to the stationary Schwarzschild 
black hole, we have observed that the decay becomes faster and the oscillation period becomes 
shorter, thus both \u)i\ and uir increase with time. 

Next, we examine the behavior of the QNMs with the increase of the charge. The Figs. 3-6 
display the frequency u>r and as a function of v for 1 = 2, evaluated at r = 5. The mass of 
the black hole is M(v) = 0.5 db 0.001i>. Both ujr and loj can be determined. Different from the 
stationary case, all curves of ujr and uj for different q increase or decrease linearly with respect to 
v. In general, as q increases from to 1, lor increases for different v. uj exhibits the behavior of 
first decrease and latter increase with the increase of q for different v, which was also observed in 
the study of stationary case [23] . 
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Figure 2: Temporal evolution of the field in the background of Vaidya metric (q=0) for / = 2, evaluated at 
r = 5. The mass of the black hole is M(v) = 0.5 ± O.OOlu. The field evolution for M(v) = 0.5 + O.OOlu and 
M(v) = 0.5 — O.OOlw are shown as the top curve and the bottom curve respectively. For comparison, the 
oscillations for M — 0.5 is given in the middle line. 
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Figure 3: The frequency ojr for 1 — 2, evaluated at r — 5 in the linearly absorbing black hole M(v) — 
0.5 + O.OOlw. 

In order to compare with the QNMs in the stationary RN black hole, we plot the curves 
CR(l,q)/uJR and ci(l,q)/uji with respect to time v as shown in Fig.7-10. cn(l,q) and ci(l,q) are 
defined in It is interesting to note that all these curves for M(v) = 0.5+0.00Tv are nearly equal 
for different q. Similar results are obtained for the curves corresponding to M(v) = 0.5 — O.OOlw. 
The slope of the curve cr(1 , q) / 'to r (or ci(l,q)/uji) is equal to the M, which suggests that the 
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Figure 4: The frequency u>i for I — 2, evaluated at r = 5 in the linearly absorbing black hole M(v) 
0.5 + O.OOlu. 
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Figure 5: The frequency ojr for I = 2, evaluated at r = 5 in the linearly evaporating black hole M(v) 
0.5-O.OOlu. 



formula 1)11(1 is still a good approximation. In fact we find 

M (« — u J — p — , M (v — u J ~ 



(29) 



ujr(v — v') u>j(v — v') 

where the parameter v' represents the retarded effect. As shown in Fig. 11, the frequencies lor for 
different position radius r with I = 2, q = and M{v) = 0.5 — O.OOlv are plotted. Clearly, the 
temporal evolution of the field \I/(r, v) at a radius r = 20 has a retard corresponding to the position 
r = 10. 
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Figure 6: The frequency u>i for 1 = 2, evaluated at r — 5 in the linearly evaporating black hole M(v) 
0.5-0.001?;. 
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Figure 7: cr(Z, q)/ujR as the function of time u for I = 2, r = 5 and M(«) = 0.5 + O.OOlw. 



Now we extend our discussion to the quick change of the black hole mass. For the linear 
model we adopted is shown in Fig. 12. With the increase of the black hole mass, we find the same 
qualitative behavior of the change of the QNM frequencies as that we observed above. The time 
scale of the variation of the black hole mass is now shorter than that of the QNM frequencies. 
The observed behavior of the quasinormal ringing is shown in Fig. 13. ujh and uji are determined 
and exhibited in Fig. 14, 15. We have shown that in the case when the variation timescale of black 
hole mass is shorter than that the QNM frequency, the waveform and the behavior of the QNM 
frequencies are the same as those we discussed for the slowly changing of the black hole mass. 

In the physical sense, the waveform and the frequencies of QNM's can be understood in terms 



12 



0.60 



0.58 



0.56 



M=0.5+0.001v 
■ q=0 

q=0.7 
a q=0.999 



. * I 



St 0.54 
o 



. *a" 



0.52 



0.50 




100 



Figure 8: c/(/, q)/ui as the function of time v for I = 2, r — 5 and M(u) = 0.5 + O.OOlv. 
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Figure 9: cr(1, q)/ujR as the function of time v for / = 2, r = 5 and M(u) = 0.5 — O.OOlw. 

of wave scattering in a given spacetime [27] and they carry significant information about the 
background. The amplitude of the wave decreases with the number of scatterings. The repeated 
scatterings will retard the influence on the temporal evolution of the wave by the time-dependent 
potential due to the quick change of the black hole mass, which leads the time scale for the change 
of the mass to be different from the time scale for the change of the frequencies. In general case, 
the repeated scatterings cause the time scale for the change of the QNM frequencies to be longer 
than the time scale for the quick change of the black hole mass. 

We have also extended our discussion to the model with exponentially decreasing mass and 
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Figure 10: ci(l,q)/uii as the function of time v for Z = 2, r = 5 and M(v) — 0.5 — O.OOlw. 
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Figure 11: lor as the function of time v for I = 2, q = 0, M(v) = 0.5 — O.OOlf, evaluated at different position 
r = 10, r = 20 and r — 30. The temporal evolution of the field ^(r, v) at radius r — 20 (or r — 30) has a 
retard to the position r = 10. 



(30) 



charge [26], 

mo v < f o 

M(u) = ^ mnexp(— at>) i>o < v < v i 

mi v > v% 

Q(v)=qM(v), 

where mo, mi, a and are constant parameters. Changing the sign before a, (|30[) can be used to 
describe the absorption process. In order to compare the different results between the exponential 
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Figure 12: Considering the black hole mass:? = 0.999, M(v) = 0.5 for v < 60, M (v) = 0.5 + 0.05(u - 60) for 
60 < v < 63 and M (v) = 0.65 for w > 60. 
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Figure 13: At r = 5, the observed behavior of the quasinormal ringing for I = 2. 

model and the linear model, we assume that there are two absorbing black holes. The mass of one 
black hole increases from M(vq) = 0.5 to M{v\) = 1 exponentially, while the other increases from 
M(vq) = 0.5 to M(v\) = 1 linearly. Fig.16 shows the cn(l,q)/ujR and ci(l,q)/wi as functions of 
time v for the linear and exponential model. The relative parameters in both models are selected 
as I = 2, r = 5, q = 0.999, vq = and v\ = 150. The two dot lines are the mass functions M{v — v') 
for models and (|2f?|) with v' = 8. In Fig. 17 we plotted cn(l,q)/MujR and ci(l,q)/Mui as a 
means to check the validity of (|29[) . The behaviors we have observed in Figs. 16, 17 also hold 
for other values of q (see Fig. 18 and Fig. 19 for q=0 as an example). In the linear model both 
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Figure 14: Change of the real part of the quasinormal frequency. Black hole mass starts to increase when 
v — 60, and the real part of the QNM frequency starts to decrease at the same moment. It is clear that time- 
scale of the change of the QNM frequency is much longer than the change of the black hole mass(c?w = 3). 
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Figure 15: Change of the imaginary part of the QNM frequency. Black hole mass starts to increase when 
v = 60, and the imaginary part of the QNM frequency starts to increase at the same moment. It is clear 
that time-scale of the change of the QNM frequency is much longer than the change of the black hole 
mass(du = 3). 

1/ujr and depend linearly on v, while in the exponential model both 1/ur and l/uoj depend 
exponentially on v. The consistency of M(v — v') and cn(l,q)/ujR (or ci(l,q)/ui) implies the 
validity of the equations (|29|) . We have also checked the validity of the equations (|29|) for different 
I, such as I = 3, 4, 5. The result keeps the same when the variation time scale of the black hole 
mass is shorter than that of the QNM frequency. 
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Figure 16: CR(l,q)/u>R and ci(l,q)/uji as the functions of time v for the linear and exponential model. In 
the linear model, A = —1/150, and in the exponential model a = — (ln2)/150. Both models describe the 
absorption process of the black hole, where the mass increases from M(vq) — 0.5 to M(v\) = 1. The relative 
parameters in both models are selected as I = 2, r = 5, q = 0.999, v = and v\ — 150. The two dot lines 
are the mass functions M(v — v') for models lj3U|) and (|2"tj)) with v' = 8. 
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Figure 17: cn{l,q)/MuiR and Ci(l,q) / Mloj as the functions of time v. The corresponding parameters are 
the same as chosen in Fig. 16. 

V. SUMMARY AND DISCUSSION 

We have studied the evolution of the massless scalar field propagating in a time-dependent 
charged Vaidya black hole background. The two generalized tortoise coordinate transformations 
were used to study the evolution of the massless scalar field. In the first case we used both the r + 
and r_, which is appropriated for the case of large charge in numerical calculation. In the other 
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Figure 18: CR(l,q)/u)R and Ci(l,q)/u>i as the functions of time v in the case of q=0. The corresponding 
parameters are the same as chosen in Fig. 16. 
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Figure 19: cn(l,q)/MtUR and ci(l,q)/Muji as the functions of time v. The corresponding parameters are 
the same as chosen in Fig. 16. 

case we only used the r + , which is appropriated for the case of small charge. The wave equation 
with a time-dependent scattering potential is derived. In our numerical study, we have found the 
modification of the QNMs due to the temporal dependence of the black hole spacetimes. In the 
absorption process, when the black hole becomes bigger, both the imaginary frequency \uji\ and 
the real frequency ojr decrease with time. However, in the evaporating process, when the black 
hole loses mass, both \ojj\ and ur increase with the increase of time. The study shows that, for the 
slowest damped QNMs, the approximate formulas in stationary RN black hole is a good description 
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for the time-dependent charged Vaidya black hole, which can be expressed as (|'29|) . 

In general, we conclude that the time dependent processes are well described by an adiabatic 
approximation with a delay in the definition of the mass function, borrowing the results at each 
moment from the problem of a black hole of a time- independent masses, instantly equal to the 
map of the equivalent Vaidya black hole at a retarded time. This implies a huge simplification for 
the problem of determination of signals from time-dependent processes near black hole. Finally we 
notice that the period corresponding to QNM is much smaller than the time of flight of a relativistic 
object in the black hole, actually by several orders of magnitude. We thus conclude that the value 
we used for the A parameter is large enough to include realistic cases. 
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